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ABSTRACT 

Context. The diagnosis of new high-resolution spectropolarimetric observations of solar prominences made in the visible and near- 
infrared mainly, requires a radiative modelling taking into account for both multi-dimensional geometry and complex atomic models. 
Aims. Hereafter we contribute to the improvement of the diagnosis based on the observation of He i multiplets, by considering 2D 
non-LTE unpolarized radiation transfer, and taking also into account the atomic fine structure of helium. 

Methods. It is an improvement and a direct application of the multi-grid Gauss-Seidel/SOR iterative scheme in 2D cartesian geometry 
developed by us. 

Results. It allows us to compute realistic emergent intensity profiles for the He i /1 10830 A and D3 multiplets, which can be directly 
compared to the simultaneous and high-resolution observations made at THeMIS. A preliminary 2D multi-thread modelling is also 
discussed. 
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. 1. Introduction 

Solar prominences (filaments) are made-up of dense and cool 
chromospheric plasma hanging in the hot and low density corona 
(Tandberg-Hanssen 1995 ). Besides its intrinsic interest as a nat- 
I ural laboratory for plasma physics, the study of these structures 
. is also of general interest in the frame of space weather stud- 
■ ies. Indeed, among other closed magnetic regions such as active 
regions, eruptive prominences are often associated with coro- 
nal mass ejections, or CMEs, that are huge plasma "bubbles" 
ejected from the solar corona and able to strongly affect Sun- 
Earth relationships, by their interactions with the terrestrial mag- 
netosphere (see e.g., Gopalswamy et al. 2006 for a recent review 
upon the various precursors of CMEs). 

Despite systematic observations made since the nineteenth 
century and decades of study, prominence formation mecha- 
nisms are still not well understood. In particular, yet no theory 
can fully explain their remarkable stability in a hotter and less 
dense medium. However, since the plasma JS is low in promi- 
nences, the magnetic field is very likely to play a major role 
in the physical scenarios which could explain prominences for- 
mation, stability and, finally, the triggering of these instabilities 
leading to CMEs. 

He I multiplets such as /110830 A in the near-infrared, and the 
Fraunhofer "yellow line" D3 at /15876 A are the best tools, so far, 
to study prominence magnetic fields. Indeed, only a few spectral 
lines are intense enough for ground-based observations in the 
optical spectrum of solar prominences i.e., at these wavelength 
at which spectropolarimetry is usually done. These helium mul- 
tiplets provide, even if they are fainter than Ha for instance, 
the most suitable information necessary for the purpose of de- 
termining the magnetic field pervading the prominence plasma, 
even though the helium spectrum, as observed at high-spectral 
resolution, reveals atomic fine structure. First spectropolarimet- 



ric observations of prominences and associated results about 
the magnetic field properties have been reviewed by Paletou & 
Aulanier ( 2003 ), for instance. 

More recently, the first full- Stokes and high- spectral reso- 
lution observations of the Hei D3 multiplet made at THeMIS 
(Paletou et al. 2001) have led to a revision of magnetic field in- 
version tools (Lopez Ariste & Casini 120021) . In particular, these 
authors demonstrated how taking into account all Stokes param- 
eters, and not only linear polarization signals, aflFects the reli- 
ability of the inversion process. Furthermore, measurements of 
the ratio between the two peaks resulting from the helium atomic 
fine structure of the He i /110830 A and D3 multiplets (see Fig. 1 1 
in Lopez Ariste & Casini 2002" for instance) are often in con- 
tradiction with the commonly used hypothesis of optically thin 
multiplets (Bommier l 19771) . 

Besides, the most recent radiative models (Labrosse & 
Gouttebroze 120011 120041) still assume mono-dimensional (ID) 
static slabs and no atomic fine structure for the He i model-atom, 
which lead to the synthesis of unrealistic gaussian profiles. It is 
therefore important to use the best numerical radiative modelling 
tools in 2D geometry, as a first step, and a more detailed He i 
atomic model in order to improve our spectral diagnosis capabil- 
ity. It is also a first application of the new 2D radiative transfer 
code we have recently developed (Paletou & Leger 2007 , Leger 
et al. l2QQ7l) . 

We describe in Sect. 2 our prominence and our atomic 
model, as well as our numerical strategy. Then, in Sect. 3 we 
compare and validate our results against previous works of 
Labrosse & Gouttebroze (2001, 200 4). We also show how geo- 
metrical eflfects can influence line profiles. In Sect. 4 we present 
2D emergent intensity profiles for the Hei /110830 A and Di, 
multiplets which can be compared to high-resolution spectropo- 
larimetric observations made at THeMIS. Finally, we discuss 
some preliminary results of 2D multi-thread radiative models. 
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2. 2D non-LTE radiative modelling 

2.1. Geometry and external Illumination 

We have adopted the prominence geometrical model of 
Vial (119821) . As illustrated in Fig. [T] it consists in an isolated 
2D slab standing above the photosphere. The freestanding slab 
is supposed to be homogeneous, static, isothermal and isobaric. 
It is characterized by two geometrical dimensions: its horizontal 
thickness Dy and its vertical extend (the third dimension is 
infinite in 2D), and its altitude above the photosphere H^. The 
prominence material is then assumed to be composed of neutral 
and ionized hydrogen, and helium. 

The slab is illuminated from below, symmetrically on its 
sides and bottom, by a photospheric and chromospheric incident 
radiation field. The latter is diluted according to the chosen al- 
titude above the solar surface, as explained by Paletou (1996). 
In addition, for the helium case, we consider a constant coro- 
nal illumination on sides and top surfaces. Indeed, as shown 
by Andretta & Jones (1997) and Mauas et al. (12005b . coro- 
nal illumination plays a very important role in the formation 
of the He I multiplets, through the so-called "photoionization- 
recombination" (PR) mechanism. For the sake of taking accu- 
rately into account for these external illumination conditions, we 
have chosen depth points logarithmically spaced away from the 
boundary surfaces towards the slab center, and symmetrically 
distributed. 



2.2. Numerical strategy 

The numerical strategy used here is similar to the one described, 
in ID, by Labrosse & Gouttebroze (2001 hereafter LGOl): 

- we solve the statistical equilibrium equations (SEE) includ- 
ing ionization balance, and the non-LTE radiative trans- 
fer equations self-consistently, for the multilevel hydrogen 
atom, using our 2D multi-grid GS/SOR iterative scheme 
(Paletou & Leger .2007. Leger et al. ^TMT). We thus ob- 
tain bound-level populations and electron densities (see also 
Paletou 1995, HeinzelT995 for the treatment of the ioniza- 
tion equilibrium). This sets physical conditions which shall 
be used as an input for the computation of the helium atom. 
The ratio between the helium and hydrogen total population 
densities is a = 0.1. We consider further that the helium 
atom is not an electron donor in the prominence, so that the 
electron density given by the hydrogen computation is kept 
constant for the helium computation. Finally the Lyman con- 
tinuum opacity and its emissivity are used to treat the helium 
ultra-violet (UV) lines pumped by the hydrogen continuum. 

- Then, we solve self-consistently the SEE and the radia- 
tive transfer equations for the multilevel helium atom, us- 
ing the same basic iterative scheme. To take into account the 
atomic fine structure, we consider the overlapping of all fine 
structure transitions into a single multi-component radiative 
transition. For example the ^3 multiplet is a combination 
of 5 fine structure transitions (see e.g.. Fig. 3 in House & 
Smartt [T982l) . We thus use, for the helium case, the Rybicki 
& Hummer (1992) full-preconditioning strategy in order to 
solve the SEE including overlapping transitions. 

In our radiative transfer codes for the hydrogen and helium 
atoms, we assumed complete redistribution in frequency (CRD). 




Fig. 1. 2D geometrical prominence model: the isolated, isother- 
mal and isobaric slab is standing horizontally above the photo- 
sphere. For hydrogen and helium atoms, we consider a chromo- 
spheric and photospheric incident radiation coming from below, 
symmetrically on its sides and bottom. For the helium case, we 
consider in addition a constant coronal illumination coming from 
above onto the slab. 



2.3. Atomic models 

We adopte d the 5 -level plus a continuum hydrogen atom model 
of Paletou (11995b , whose data is consistent with earlier models 
of Gouttebroze et al. flWT). 

For neutral helium, we considered two diff'erent cases: 19 
bound-levels, up to ^ = 4, plus a continuum, with no atomic fine 
structure (hereafter our Hen4 model), or 17 bound-levels, up to 
n = 3, plus a continuum, with the atomic fine structure of 2^P, 
3^P and 3^D levels (hereafter, our Hen3sf model). 

Energy levels and statistical weights were taken from the 
NIST database (Ralchenko et al. |2Q^ . All other atomic data 
are consistent with those of LGOl: eff"ective collision strengths, 
collisional ionization coefficients and spontaneous emission co- 
efficients are from Benjamin et al. ('1999'), collision strengths not 
defined there are from Benson & Kulander ( 1972 ), and photoion- 
ization cross sections are from TOPbase (Fernley et al |1987l ). 

Our radiative transfer code can explicitely consider colli- 
sional rates among them and to/from all fine structure sublevels. 
For our Hen3sf model, we set to the rates between sublevels of 
the same level. Otherwise, we used, for each sub-level, the rates 
given by the later authors even though they were summed over 
the fine structure. We are aware of this inconsistency which may 
have only a very small influence on the results exposed hereafter. 

Chromospheric illumination for all radiative transitions were 
taken from Heasley et al. ( 1974) , with some additional details 
which can be found in Labrosse et al. (120071) . Coronal illumina- 
tion for UV lines and continua were taken from Tobiska ([1991]) 
and Walhstr0m & Carlsson (.1994) . 



3. Preliminary results 

As a first step, results obtained with our 2D radiative transfer 
code for the hydrogen atom are in very good agreement with 
those of He asley & Milkey (1983), Gouttebroze et al. (1993D 
and Paletou (IT995lP l. 

As a second step, we have compared the results obtained 
with our upgraded 2D radiative transfer code for the helium 



^ In this article, partial redistribution in frequency (PRD) was used 
by default, but we could make direct comparisons with a CRD version 
of this antique code. 
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Fig. 2. Mean population and electron densities as a function of 
temperature for two pressures, pg = 0.02 (top) and 0.2 (bottom) 
dyn cm"^. Numbers in the upper-right frames refer to the energy 
level as follows: 1 for the ground state of He i and 20 the ground 
state of He ii. The mean population densities for these two levels, 
as well as for the electron density, are divided by 10^. The other 
codes for levels are: 4 for 2^P, 5 for 2^P,9 for i^D and 10 for 



atom with those of LGOl, as these last authors have made com- 
prehensive comparisons with the pionneering works of Heasley 
et al. (1974) and Heasley & Milkey (1978). In order to do so, 
we have used our Hen4 model, without atomic fine structure. 
Also, their radiative code is ID plane parallel so that they have 
neither coronal illumination on the top surface, nor dilution ef- 
fects varying with the altitude, as in 2D. However, their helium 
atomic model is more detailed than ours, with 29 bound-levels 
up to ^ = 5 for He i, 4 bound-levels for He ii, and He iii. Unlike 
us, they also take into account PRD for the resonance lines H i 
\.ya and LyyS, He i ^584 A and He ii .1304 A. 

We have modeled a 2D prominence with = 100 000 km 
and Dy = I 000 km. Its temperature, T, could take the respec- 
tive values: 6 000, 10 000, 14 000, 17 000 and 20 000 K, and the 
gas pressure pg could take 2 values: 0.02 and 0.2 dyn cm"^. The 
micro turbulent velocity was fixed at ^ = 5kms"^ The bottom 
of the slab was set at Hq = 10000 km above the solar surface. 
We have chosen a 2D spatial grid with 123 x 123 logarithmically 
spaced points. We have also used Doppler profiles monotoni- 
cally sampled with a 0.1 step in Doppler width units. 

In Fig. [21 mean population densities variations against tem- 
perature are displayed. Mean densities are defined here as: 




- 14000 
■ 17000 
20000 



.0 0.8 0.6 0.4 0.2 0.0 0.2 0.4 0.6 0.8 

ae+4 10830 A 



- 6000 

- 10000 




-1.5 ^1.0 -0.5 0.0 0.5 1.0 1.5 
AX (A) 



Fig. 3. Emergent intensities for lines Hei A5S4 A, /15876 A, 
/16678 A and /110830 A, for a pressure pg = 0.02 dyn cm'^, 
for 5 diff'erent temperatures, and for = 100000 km and 
Dy = I 000 km. These intensities are computed at slab mid- 
height, z = for a line-of-sight perpendicular to the slab. 




Ni = 



X""^ ni(y,DJ2)dy 



A; 



Fig. 4. Same as Fig.[3]for a vertical geometrical extension = 
10000 km. 



where nfy, 0^/2) is level / population at the position y at mid- 
height of the slab (z = 0^/2). We have focused on the popula- 
tions of 5 energy levels of Hei: 1^^, 2^P, 2^P, 3^D, 3^D labelled 
1, 4, 5, 9 et 10, and the ground state of Hen labelled 20. For the 
sake of clarity, mean population densities of the ground states of 
He I and He ii, as well as the electron density, are divided by 10^ 
on the figure. 

In the low pressure regime, pg = 0.02 dyn cm"^, we recover 
the same general behaviour as the one shown by in LGOl: an 
increase in temperature reduces the neutral helium mean popu- 
lation and increases the ionized helium mean population, until 
the latter becomes greater than the former. This transition tem- 
perature is around 8 000 K for our model, whereas it is 11 000 K 
in LGOl. This is explained by PRD eff'ects, as LGOl models in 
CRD give the same result as ours (N. Labrosse, private commu- 
nication). 

For larger pressures, pg = 0.2 dyn cm"^, it is shown in Fig. [2] 
that all mean population densities increase with temperature but 
the ground state of He i. As emphasized indeed by LGOl, the he- 
lium ionization is less important than in the low pressure case. 
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because of the optical thickness of the helium continuum which 
prevents incident radiation from reaching the core of the slab. 
We also recover this effect using our 2D code, when the verti- 
cal extension of the slab is made large enough, mimicking a ID 
vertical slab model though. 

We could also check a posteriori, comparing n(Heii) and ne, 
that the assumption that He is not an electron donor remains 
valid for most of the range of parameters we used. It becomes 
however questionable for both high pressure and high tempera- 
ture cases, as shown in Fig.O 

In Fig. [3l emergent intensities for lines Hei A5S4 A, 
/15876 A, ^6678 A and .110830 A are displayed for a vertical 
extension = 100000 km, for a pressure pg = 0.02dyncm"^ 
and for five diff'erent temperatures. Intensities are computed at 
mid-height of the slab, at z = 0^/2, and for a line-of-sight (los) 
perpendicular to it. 

The general behaviour of these intensities with temperature 
is in excellent agreement with results presented by LGOl (see 
e.g., their Fig. 8). Our values for the intensity of the resonance 
line He I ASM A are smaller than those of LGOl, whereas they 
are greater for the optically thin /15876 A and /16678 A lines. 
This could be explained by PRD eff'ects which play an important 
role in determining the shape of the resonance line He i /1584 A 
as well as in the He i levels 1 and 4 populations. These PRD ef- 
fects have consequently an impact on the line Hei /15876 A as 
its lower level is the level 4, and on ionization via the neutral 
helium continuum A504 A. The ground state of He ii mean pop- 
ulation density is smaller in LGOl. As the triplet levels are pop- 
ulated through the PR mechanism (see Andretta & Jones 119971 
and references therein) from the ground state of He ii, our D3 
multiplet emergent intensity is greater than the ones of LGOl. 

In Fig. m the same emergent intensities are displayed for a 
smaller vertical extension of Z)^ = 10000 km (note also that 
we used, in both case, the same fixed dilution factor for the inci- 
dent radiation, in order to put in evidence the only eff'ects due to 
diff'erent vertical slab extensions). Emergent intensity values are 
larger, with relative diff'erences ranging between 15 and 40%, 
than those displayed in Fig. [3] This shows that significant 2D 
geometrical eff'ects take place when one reduces the geometri- 
cal extension of the slab, and how they impact the magnitude of 
emergent intensities. Such kind of eff'ects were first put in evi- 
dence on Ha by Paletou (.1997) . 

4. A realistic spectral synthesis of Hei multiplets 

We now focus on Hei D3 and /110830 A multiplets, as our pri- 
mary aim is to model realistic emergent line profiles consistent 
with high-resolution spectropolarimetric observations made at 
the solar telescope THeMIS (see e.g., Paletou et al. 120011) . With 
such an instrument, the two "red" and "blue" peaks are, indeed, 
clearly resolved (see also e.g., Merenda et al. 2006, Trujillo 
Bueno et al. 2002, for data taken at the German VTT with the 
TIP polarimeter, in the near-infrared). 

Hereafter, we shall use in our 2D radiative transfer code the 
Hen3sf atomic model which takes into account the atomic fine 
structure for the 2^P,3^P and 3^D levels of He i. 

4.1. Single slab models 

We have modeled a 2D isobaric, isothermic and static promi- 
nence with a vertical extension = 30000 km and a hori- 
zontal extension Dy = 5 000 km. Its temperature is T = 8 000 
K, and the gas pressure is pg = 0.05 dyn cm"^ which are rather 
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Fig. 5. Normalized emergent intensities for the Hei /110830 A 
multiplet obtained with the Hen4 atomic model (dashed lines), 
with the Hen3sf atomic model (thick line), and from THeMIS 
observations of June 2007, vs. the wavelength shift with re- 
spect to the multiplet central (vacuum) wavelength defined as 
A = 10832.7 A. Synthetic profiles are for a line-of-sight per- 
pendicular to the slab. Both the fine- structured and the observed 
profiles were normalized to the maximum amplitude of the Hen4 
synthetic profile. The observed one was also slightly shifted in 
frequency, for the sake of comparison with the Hen3sf one. 
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Fig. 6. Same as Fig. [5] for the Hei D3 multiplet whose central 
(vacuum) wavelength is 5877.4 A. 



typical values for the modelling of quiescent prominences (see 
e.g., Gouttebroze et al. 1993). The microturbulent velocity is 
^ = 5 kms"^. The bottom of the slab is set at Hq = 10000 km 
above the solar surface. We have chosen a 2D spatial grid with 
243 X 243 logarithmically spaced points. 

In Figs. [5] and [6l emergent intensities for Hei /110830 A 
and D3 multiplets are displayed. They were computed at the 
slab mid-height, z = 0^/2, for a los perpendicular to the slab, 
and the specific intensity was normalized to the maximum value 
obtained with model Hen4. Then, maximum intensity values 
are 21871 erg s"^ cm'^ sr"^ A"^ for the /110830 A multiplet and 
6473 erg s"^ cm"^ sr"^ A"^ for the D3 multiplet. Taking into ac- 
count the atomic fine structure, we obtain indeed, 2D emergent 
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Fig. 7. A sketch of the 2D multi-thread prominence model we 
have adopted. Individual threads are superimposed along the 
line-of-sight, and there is no radiative interaction between them. 
Emergent intensities are computed as formal solutions of the ra- 
diative transfer equation along a line-of-sight perpendicular to 
the threads vertical extension, in order to take into account the 
contribution of each individual thread. 



line profiles which are directly comparable to high-resolution 
spectroscopic observations, characterized by two well-resolved 
subcomponents. 

We define Ir as the ratio between the intensity of the largest 
subcomponent peak and the intensity of the smaller peak, that 
is, hiuelhed for the Z)3 multiplet and Iredlhiue for the /110830 A 
multiplet. We remind here that, when the He i multiplet be- 
comes optically thick, the Ir ratio is lower than 8 (House & 
Smartt 1982, Landi Degl'Innocenti 1982). Furthermore, statis- 
tics made on this ratio give a mean value of 6 (Lopez Ariste & 
Casini 2002), which is not consistent with the commonly used 
hypothesis that the multiplet is optically thin. 

We can evaluate the optical thickness of the two multiplets 
obtained with our 2D model. At the slab mid-height, z = 0^/2, 
and for the frequency corresponding to the maximum intensity 
value: 

t(10830) = 5, 1 X 10"^ 
t(Z)3) = 4,4 X 10"^ 

The two multiplets are thus optically thin, and we also find 
that t(10830) ~ 10 x t(Z)3), as pointed out by Andretta & 
Jones d 19971) . In these two cases, Ir ratios are indeed around 8. 

Having performed several tests with temperatures varying 
from 6000 to 20000 K and for gas pressures ranging from 
0.02 to 0.2dyncm"^, we found that the two multiplets are op- 
tically thicker for both high pressure and high temperature at- 
mospheres. However, even in these regimes, the optical thick- 
nesses of the D3 multiplet never go beyond r ~ 0.1. These 
results are also in accordance with ID models of Labrosse & 
Gouttebroze (120041) . Eor instance, these authors found no model 
for which the D3 multiplet is optically thick, and t( 10830) > 1 
for models where T > 20000 K and pg > 0.64dyncm"^ i.e., 
conditions which are not representative of those expected in qui- 
escent prominences. 

Which kind of model could reproduce the observed ratio be- 
tween subcomponents of the 557.6 and 1083. nm multiplets of 
He I? 

4.2. 2D multi-thread models 

We have thus decided to model prominences diff'erently, using 
the fact that they are made of small-scale threads. This is sup- 
ported both by observations (e.g., Berger et al. 120081 Lin et 
al. I2005>) and by MHD simulations (e.g.. Low & Petrie 12005 ). 
The effect of an increased penetration of the incident radia- 
tion, using multiple layer modelling, was also shown in ID by 
Gouttebroze et al. ([2Q02D. Also, such a 2D non-LTE spatial fine- 
structure radiative modelling of solar prominences was recently 
used by Gunar et al. (2007), for the purpose of interpreting UV 
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Fig. 8. The Ir ratio (solid line) and the optical thickness inte- 
grated on the slab height (dashed line), for the Hei /110830 A 
multiplet versus the number of individual threads, for a gas pres- 
sure Pg = 0.5 dyncm"^ and a temperature 7 = 8 000 K. 




number of threads 



Fig. 9. Same as Fig.[8]for the Ir ratio and the optical thickness of 
the He i multiplet. 



observations of the Lyman series of hydrogen made with the 
SUMER spectrograph on-board SoHO. 

Hereafter, all the 2D threads (see Fig. [7]) are supposed to be 
identical and located at the same altitude, with similar horizon- 
tal and vertical extensions Dy = = I 200 km which cor- 
responds to a 1.5 arcsec angular resolution compatible with 
what can be decently achieved while doing spectropolarimetry 
of prominences using ground-based observatories. Such a choice 
can obvisouly be debated, and we could have considered (much) 
smaller threads, beyond spatial resolution (see e.g.. Vial 2 0061 
Heinzel 2007 ). However, we leave more detailed investigations 
about such an issue to further studies. 

We considered also that there is no radiative interaction be- 
tween individual threads. Thus, as a first step, we applied our 2D 
radiative transfer code for one individual thread. Then we com- 
puted emergent intensities as the formal solution of the radiative 
transfer equation along a los perpendicular to each thread's ver- 
tical extension, in order to take into account the contribution of 
a bunch of threads. Therefore, the total optical thickness along 
the los is simply the sum of the optical thicknesses of each indi- 
vidual thread. 
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For each (identical) thread, we adopted a gas pressure p = 
0.5 dyncm"^ and a temperature 7 = 8 000 K corresponding to 
a commonly used temperature for quiescen t prom inences (see 
e.g., Engvold et al. 1 19901 Tandberg-Hanssen I 1995b . 

In Figs. [8] and [21 we have drawn the Ir ratio versus the num- 
ber of threads along the los, for the He i /110830 A and mul- 
tiples. For the calculation of this ratio, we considered the emer- 
gent specific intensity integrated over as 



i(y) 



/(y, z)dz 



The total optical thickness, integrated the same way, is also dis- 
played on the same plots. 

We can thus determine the number of threads required, with 
such models, to obtain a Ir ratio around 6 for the He i D3 mul- 
tiples In that case, a good value is around 30 threads, which 
corresponds, however, to a very large total geometrical width of 
36 000 km for the prominence. The corresponding optical thick- 
ness is t(Z)3) ~ 0.2. For a temperature r = 17 000 K, corre- 
sponding to a prominence to corona transition region (PCTR, see 
e.g., Fontenla et al. 1996, Anzer & Heinzel 1999 and references 
therein), the number of threads falls to 15, and riDs) ~ 0.7. 

Even though it is still preliminary, yet we favour multi-thread 
models in order to explain a number of properties of observed 
He I spectral lines. 

However, we are also aware of the fact that, even in a non- 
magnetic field regime. Stokes / is aff'ected by atomic polarization 
(R. Casini, private communication). Therefore, a more detailed 
comparison with spectropolarimetric data would require a more 
complex treatment of the statistical equilibrium. 

5. Conclusion 

We have performed here several tests in order to compare 
our new 2D radiative model against previous works. Results 
for the hydrogen atom have been compared to Heasley & 
Milkey (1983), Gouttebroze et al. (1993) and Paletou (1995). 
Results for the helium atom have been mainly compared to 
Labrosse & Gouttebroze (2001, 2004). An excellent agreement 
is reached using our 2D code in the limit of a ID geometry, with 
similar atomic models and incident radiation. 

We have shown that, taking into account the atomic fine 
structure for the 2^P, 3^P and 3^D levels of Hei in our 2D ra- 
diative transfer code, allows to compute emergent line profiles 
which can be directly compared to high- spectral resolution ob- 
servations. However, using a classical model of isothermal, iso- 
baric, homogeneous and static slabs leads to ratios between the 
subcomponents of Hei /110830 A and D3 multiplets which do 
not always support observed values, except for high temperature 
and high pressure conditions, which are not typical of quiescent 
prominences (Engvold et al. 1 1990b . 

This led us to undertake a preliminary investigation of 2D 
multi-thread modelling. It has indeed given us good hints about 
how it should be possible to reproduce observed characteristics 
of the D3 and /110830 A multiplets of Hei, such as the ratio 
between the intensity of their two subcomponents. 

Our 2D radiative model will now be more intensively used, 
taking into account smaller threads, eventually beyond spatial 
resolution (see Heinzel [20071 for a recent review) and the radia- 
tive interaction between them (e.g., Heinzel 1989). We should 
also include a PCTR for each individual thread in our model, 
since its influence on hydrogen lines have already been evalu- 
ated. 



We finally anticipate that our results will be very valuable for 
the analysis of data coming from a variety of spaceborne (e.g., 
SoHO, Hinode) and groud-based telescopes such as THeMIS. 
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